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We present correlation results for the radio halo power in galaxy clusters with the integrated thermal Sunyaev-Zel'dovich 
(SZ) effect signal, including new results obtained at sub-GHz frequencies. The radio data is compiled from several pub- 
lished works, and the SZ measurements are taken from the Planck ESZ cluster catalog. The tight correlation between 
the radio halo power and the SZ effect demonstrates a clear correspondence between the thermal and non-thermal elec- 
tron populations in the intra-cluster medium, as already have been shown in X-ray based studies. The radio power varies 
roughly as the square of the global SZ signal, but when the SZ signal is scaled to within the radio halo radius the correlation 
becomes approximately linear, with reduced intrinsic scatter. We do not find any strong indication of a bi-modal division 
in the radio halo cluster population, as has been reported in the literature, which suggests that such duality could be an 
artifact of X-ray selection. We compare the Ksz dependence of radio halos with simplified predictions from theoretical 
models, and discuss some implications and shortcomings of the present work. 
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■ 1 Introduction 

■ Precision use of galaxy clusters for measuring the cosmo- 
. logical parameters largely depends on our ability to accu- 

■ rately describe the properties of the intra-cluster medium 
, (ICM). The ICM contains about 90% of the baryonic mass 

■ in massive clusters, and reaches very high temperatures (~2- 
, 15 keV) after attaining thermal equilibrium in the deep clus- 

■ ter potential well. This thermal plasma emits in the X-rays 
, through the bremsstrahlung emission, and is also visible 

■ in the radio/sub-millimeter wavebands through the thermal 
, Sunyaev-Zel'dovich (SZ) effect: the two methods that have 

■ been highly successful in finding galaxy clusters from large- 
, area surveys and studying their detailed properties. 

In addition to this dominant thermal component the 
. ICM also hosts a rich variety of non-thermal energy sources 
and magnetic fields. The most spectacular evidence for this 
cluster-wide non-thermal phenomena comes from the ob- 
servation of giant radio halos, which are diffuse synchrotron 
emission extending over ~ 1 Mpc scales. They are pre- 
sumed rare, typically found in dynamically active systems 
(mergers), and have a smooth distribution similar to the 
ICM that indicates a connection between their powering 
mechanism and the cluster mass (e.g. Liang et al. 2000, Fer- 
rari et al. 2008). What remains unknown, however, is the 
detailed mechanism by which these GeV electrons get dis- 
tributed in the cluster volume, and whether their presence or 
absence mark special sub-populations of clusters not suit- 
able for precision cosmology. We also do not know if the 
radio halo power can be used as an accurate and unbiased 
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mass indicator in merging systems, where the standard mass 
estimation method based on hydrostatic equilibrium is ob- 
viously not adequate. 

With the goal of establishing a cleaner mass selection in 
the radio halo studies, we obtained the first radio-SZ corre- 
lation in galaxy clusters (Basu 2012). The radio data were 
collected from the literature, and for the SZ measurements 
we used the Planck Early SZ cluster catalog (Planck collab- 
oration 201 1). While we cannot answer questions regarding 
the rate of occurrence of radio halos in a given mass bin 
since we lack radio observation of a mass-limited sample, 
our work aims to provide answers to the mass scaling of 
radio halos and the scatter in their radio power. In particu- 
lar, we have attempted to find out if the observed bi-modal 
distribution between radio halo and radio "quiet" clusters, 
as seen from X-ray selected samples (e.g. Brunetti et al. 
2007), is an inherent property of radio halos, or merely a 
selection artifact. This bi-modality in the radio halo cluster 
population has been used as an argument in support of ra- 
dio halo origin models based on turbulent re-acceleration of 
electrons following cluster mergers (Brunetti et al. 2007 and 
references therein). 

Compared to X-ray luminosities, the integrated SZ sig- 
nal is a more robust indicator of cluster mass, having little 
dependence on the cluster dynamical state or other details 
of cluster physics (e.g. Motl et al. 2005, Reid & Spergel 
2006). During cluster mergers both Lx and Ksz go through 
transient variabilities, but the changes in Ys,z are less ex- 
treme than in Lx (Poole et al. 2007, Wik et al. 2008), thus 
less affecting the mass selection. The gradual late-time in- 
crease in the SZ signal should also correspond better with 
the radio halo timescales (~ 1 Gyr, derived from the spatial 
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Fig. 1 Left: P1.4 - Fsz correlation for a sample of clusters with radio halos as compiled by Cassano et al. (2007). The 
integrated SZ signal from the Planck catalog have been scaled to inside the halo radius assuming the "universal" GNFW 
pressure profile. The slope of the scaling is slightly super-linear, with intrinsic radio scatter roughly 0.28 dex. Right: 
Pi A - Ysz correlation and testing of a bi-modal division in the radio halo population using the Brunetti et al. (2009) 
sample. Here Ysz represents the total SZ signal inside radius SR^qo (as given in the Planck catalog), and the slope is steeper 
compared to the scaling inside radio halos. The red circles denote the radio halo detections, blue arrows the non-detection 
upper limits (see Brunetti et al. 2007) and the orange squares is a mini halo (A2390). The short dashed line represents the 
fit for halos only, and the long dashed line when upper limits are included. The shaded regions mark the 95% confidence 
intervals for the corresponding fits. 



extent of the radio halos). As a result, we can expect the SZ 
selection to reveal the presence of radio haloes in late merg- 
ers and other massive clusters which are potentially left out 
in X-ray luminosity limited samples. 

2 Results for the radio-SZ correlation 

In the following we present some highlights of our work 
using the Planck ESZ clusters that have known radio 
halo data in the literature. All results are derived us- 
ing the ACDM concordance cosmological parameters. The 
quantity Fsz denotes the intrinsic Compton K-parameter, 
Ydj^E{zy'^^ , where cIa is the angular distance and E{z) is 
the ratio of the Hubble parameter at redshift z to its present 
value. The Planck catalog provides the integrated Y parame- 
ter within the radius SR^oo, obtained from a multi-frequency 
matched filtering (see Planck collaboration 2011). This is 
much larger than the typical radio halo extension, so we 
also obtained a scaled Ysz value inside individual radio halo 
radius using the universal pressure profile of Arnaud et al. 
(2010). Both these total and scaled integrated SZ signals are 
then correlated with the radio halo power. 

2.1 Ysz and mass scaling of radio halos 

For radio catalogs we used the compilations from Giovan- 
nini et al. (2009) and Cassano et al. (2007), which give ra- 
dio halo power measured at 1 .4 GHz. The latter sample is 
used due to its better definition of radio halo radius. The cor- 
relation of the radio power with the total SZ signal (isz(< 
5^500), as measured by Planck ) follows a roughly quadratic 



power- law: Py oc 7^2^*'"^ for the Giovannini et al. (2009) 
sample, and Py oc Kg^'^'^^ for the Cassano et al. (2007) 
sample. Scaling the SZ signal inside halo radius makes the 
correlation approximately linear: e.g. for the Cassano et al. 
(2007) sample we obtain Py oc Yszi< ./?//)'-^*° ", with intrin- 
sic scatter 0.24 dex (FiglT]/e/f). We use this slightly super- 
linear relation for obtaining the mass scaling for radio halos, 
although we emphasize from the current data a linear rela- 
tion between radio and SZ power inside halos is a fully valid 
result. This linear relation is expected to become sub-linear 
at lower frequencies for curved spectrum of radio halos. 

The relation between the integrated SZ signal and clus- 
ter mass is Ysz °^ fganM^J^^Eiz)'^^, which has been found 
to be observationally robust (e.g. Andersson et al. 2011, 
Planck collaboration 2012). /gas is the gas mass fraction 
which varies weakly with the cluster mass in the high-mass 
regime: /gas K M°;/o'* (e.g. Bonamente et al. 2008). Us- 
ing the same gas fraction at all radii we therefore obtain 
Py oc M^g^'^'^''^ for the total virial mass, or Py oc M^'=^"'^^ for 
the mass enclosed within a radio halo. 

The mass scaling of radio halos within the halo radius 
is consistent with previous estimates using X-ray hydro- 
static masses (e.g. Cassano et al. 2007). The global scal- 
ing with M500 or Mvir is steeper than that determined from 
X-ray luminosities (e.g. Cassano, Brunetti & Setti 2006), 
and can be considered more robust due to Ysz being a bet- 
ter mass proxy. Directly applying the Lx - M500 scaling to 
our samples, we obtain Py oc M^^^y which is consistent with 
the mass scaling inside radio halos, but much flatter than 
the scaling with total virial mass. This is likely due to the 
very peaked nature of the X-ray luminosity profile, causing 
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most of the X-ray emission originating from within ~ r25oo 
(which is similar to the radio halo radius), and therefore cor- 
relating more directly with the physical state of the ICM 
inside radio halo radius. 



2.2 The issue of bi-modality 

For testing whether there are two distinct populations of 
clusters, namely those hosting powerful radio halos and 
those without, we use radio halo catalogs that include non- 
detections. A systematic study of radio halo non-detections 
in an X-ray selected sample is done by Venturi et al. (2008), 
using GMRT observation at 610 GHz. However, this sam- 
ple contains too few clusters which have Planck SZ mea- 
surements to obtain any robust fit {Fig^ left). We therefore 
use the compilation by Brunetti et al. (2009), which lists 
the same GMRT results scaled to 1.4 GHz with a list of 
other halo detections. In their work, non-detection upper 
Umits were obtained by simulating fake radio halos in the 
GMRT data and scaling to 1.4 GHz by using a = 1.3 for 
the spectral index. The result of the correlation is shown in 
FigdJr/g/ir. There appears to be a hint of bi-modal division, 
but the number of radio halo non-detection clusters in the 
Planck catalog is too low to confirm this with certainty. The 
point to note here is, we do not see high-Ksz objects with 
radio halo non-detections, as in the case of highly X-ray lu- 
minous radio "quiet" cool core clusters. A similar result is 
obtained when the radio halo power is correlated against the 
total gas mass inside R500, by dividing the SZ signal with the 
mean temperature: Mgas = i'sz(^50o)/7'x HempHcrT jnieC-). 
We therefore conclude that the strong bi-modal division ob- 
served in the X-ray selected cluster samples is likely due to 
the inclusion of X-ray bright cool core clusters with lower 
masses, which do not show extended radio halos. 

An interesting comparison between the X-ray and SZ 
selections can be obtained if we compare the relative rate 
of occurrence of radio halo clusters in the original X-ray 
selected samples with that in the Planck catalog (Table 
1). Both the samples of Venturi et al. (2008) and Rudnick 
& Lemmerman (2009), which are complete above certain 
X-ray luminosity thresholds, have radio halo fraction of 
roughly 20%. This fraction increases to roughly 50% when 
we look for the counterparts in the Planck catalog for these 
samples. We emphasize that this is not an objective test of 
radio halo statistics from SZ selection, since none of these 
samples are selected a priori based on their SZ signal, hence 
they are not complete above a mass threshold. Neverthe- 
less, the point to highlight here is that the fraction of radio 
halo clusters is clearly increasing in an SZ-based sample, 
as opposed an X-ray luminosity limited samples. The most 
likely explanation is that Ksz is a better mass proxy and the 
Planck ESZ catalog is specifically showing the most mas- 
sive clusters in the universe which are likely to harbor a sig- 
nificant non-thermal electron population. 
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Table 1 Relative frequency of finding radio halos in dif- 
ferent X-ray selected cluster catalogs, and the correspond- 
ing rate of occurrence in the Planck ESZ catalog. 



Sample 


no. of 


of which 


Planck 


of which 




clusters 


radio halos 


detects 


radio halos 


V08 


26 


6 


9 


5 


B09* 


44 


21 


20 


16 


R09 


72 


14 


27 


12 



Samples: V08=Venturi et al. 2008; B09=Brunetti et al. 2009; 
R09=Rudnick & Lemmerman 2009. * Not X-ray complete 



2.3 Correlation at low frequencies 

The question of radio halo statistics and their mass scaling 
at sub-GHz frequencies is highly relevant for the ongoing 
and future radio surveys like LOFAR and SKA. We tried to 
obtain similar radio-SZ correlations for low-frequency ra- 
dio halo observations of X-ray complete samples, but the 
current data sets do not allow robust constraints. Venturi et 
al. (2008) sample is based on the GMRT observation at 610 
MHz of a X-ray selected sample. Out of 26 clusters only 9 
are present in the Planck catalog, which does not allow a 
robust regression analysis, but the result highlights the scat- 
ter in the radio halo power and inter-mixing of halos and 
non-detections (Fig|2]/e/f). We remark, however, that inad- 
equate sampling of the short baselines in the GMRT 610 
MHz measurements can result in serious under-estimation 
of radio halo fluxes (e.g. Dallacasa et al. 2009). 

Rudnick & Lemmerman (2009) re-analyzed the 
WENSS survey data at 327 MHz (Rengelink et al. 1997) for 
an X-ray selected sample. The shallowness of the WENSS 
data, which is roughly an order of magnitude less sensi- 
tive than typical VLA radio halo observations scaled to this 
frequency, does not allow a comprehensive test for the bi- 
modality. However, the method used by these authors to ex- 
tract radio halo fluxes, based on simulating fake halos in 
empty fields, safeguards against the flux loss issue. Cor- 
relation of the 327 MHz radio halo power with the global 
Fsz signal (Fig^right) yields a slope much shallower than 
those observed in the 1.4 GHz samples when only > 3cr 
halo detection are considered: P0.33 °^ Y^:^^"^^-^^. The large 
statistical errors resulting from the shallow WENSS data 
hides the intrinsic scatter in the radio power, which domi- 
nates the 1.4 GHz radio uncertainties. When non-detections 
are included with the measured values, the slope becomes 
compatible with the 1 .4 GHz result at roughly 1 .2cr, but the 
scatter also increases to almost twice as much. 

The turbulent re-acceleration model of radio halo ori- 
gin predicts a synchrotron power-law spectrum that gradu- 
ally steepens at higher frequencies (e.g. Brunetti & Lazar- 
ian 201 1). Based on such steepening one can naively expect 
the radio-SZ correlation to become flatter at lower frequen- 
cies, as the radio under-luminous objects at 1.4 GHz will 
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Fig. 2 Left: Radio-SZ correlation result at 610 MHz from the cluster sample observed by Venturi et al. (2008). This 
study is based on 26 clusters selected above Lx > 5 x 10"*"^ erg/s, but the corresponding number of counterparts in the 
Planck catalog is too few to obtain any robust SZ-radio correlation. This plot highlights the intrinsic scatter in the radio halo 
power, although flux under-estimation in the GMRT data can be an issue. Right: Radio-SZ correlation at 327 MHz from the 
sample of Rudnick & Lemmerman (2009), based on WENSS radio data. The survey data is shallow, and correspondingly 
the non-detection upper limits are not sufficiently below the detection line to test a bi-modality. All data points and lines 
have the same meaning as in FiglT] 



have a larger relative gain in their flux; although Cassano 
(2010) predicts a steepening of the Py - Lx correlation at 
low frequencies based on a model population of ultra-steep 
spectrum radio halos in low-mass clusters. The correlation 
result from Rudnick & Lemmerman (2009) shows no clear 
evidence of steepening, and is at best consistent with our 
1 .4 GHz results if we consider the halos and no-detections 
as one population. The gradually steepening spectrum of ra- 
dio halos should also lead to a reduced intrinsic scatter in the 
radio power at low frequencies, which unfortunately can not 
be tested from these data sets. 

3 Discussion 

The observed linear correlation between the radio and SZ 
power inside radio halos can be provisionally explained us- 
ing simplified hadronic models, which predicts Py/isz 
Xqr n cTpp /b (e.g. Kushnir et al. 2009). Here Xcr is the ra- 
tio between cosmic ray pressure and thermal pressure, n is 
the gas density, cTpp is the p-p collision cross-section, and 
fs is the volume filling factor for magnetic fields. This sim- 
ple formulation requires constant cosmic ray proton den- 
sity and large magnetic field strength over the entire ra- 
dio halo volume, which may be unrealistic. The turbu- 
lent re-acceleration of particles for creating radio halos re- 
quires a complex mechanism, and again we use a simplified 
model based on the works by Cassano & Brunetti (2005) 
and Cassano et al. (2007). these authors obtain a scaling 

1 

Py oc YhT^J", which is in good agreement with the ob- 
served mild super-linear correlation. However, the analytic 
relationship between the radio halo power and the cluster 
mass and temperature in this model is based on several sim- 
plifying assumptions, and the parameters require fine tuning 



to match observations. The true explanation of the observed 
quasi-linear correlation between the radio and SZ signals 
can therefore be beyond these simplified models. 

The small number of non-detections in the Planck ESZ 
catalog does not allow us to conclude whether the bi- 
modality seen in the X-ray selected samples is absent - or 
merely weaker - when measured against the cluster SZ sig- 
nal. Likewise, we have not been able to address the true 
fraction of radio halos in clusters in a given mass range, 
since the current work is based on prior X-ray selected clus- 
ters that are present in the Planck catalog. Some of these 
questions will be addressed in a future publication utilizing 
public radio survey data and complete cluster samples in the 
SZ and X-ray (M. Sommer & K. Basu, in preparation). 
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